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This  paper  focuses  on  the  sizing  optimization  of  a  thermoelectric  generator,  combined  with  a  potential 
heatsink,  that  harvests  heat  from  the  human  body.  To  ensure  success,  energy  models  describing  a  ther¬ 
moelectric  module  and  a  heatsink  are  introduced;  these  two  models  are  then  simplified  for  small  temper¬ 
ature  differences  while  respecting  the  continuity  between  them  (both  with  and  without  the  heatsink).  In 
order  to  maximize  the  power  produced  by  such  a  system,  the  dual  thermal  and  electrical  matching  of  the 
generator  with  its  environment  is  studied  in  detail.  Lastly,  to  optimize  the  sizing,  a  set  of  optimization 
parameters  are  defined,  namely  leg  length,  the  capture  surface  area  of  the  thermoelectric  generator 
and,  if  a  heatsink  is  present,  its  height.  The  results  derived  help  draw  a  conclusion  on  the  pertinence 
of  adding  a  heatsink,  for  the  purpose  of  maximizing  electrical  power  while  minimizing  total  volume. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

As  part  of  a  study  on  multisource  generator-harvesting  re¬ 
sources  within  the  human  environment  [1]  (heat,  light,  move¬ 
ment),  this  paper  focuses  on  the  potential  for  thermoelectric 
generation  from  human  body  heat.  The  goal  herein  is  to  maximize 
thermoelectric  productivity  for  the  purpose  of  charging  an  accu¬ 
mulator  that  is  capable  of  supplying  a  consumer  electronic  device, 
e.g.  communicating  sensor. 

To  this  end,  the  thermoelectric  generation  [2,3]  consists  of  plac¬ 
ing  one  side  of  a  thermoelectric  module  in  contact  with  human 
skin  and  the  other  in  direct  contact  with  ambient  air,  possibly  via 
a  heatsink.  The  temperature  difference  between  the  two  sides  is 
then  used  to  directly  convert  heat  into  electricity;  this  difference 
however  remains  relatively  low,  mainly  as  a  result  of  the  small 
temperature  differential  between  the  two  heat  sources  (human 
body  and  ambient  air),  as  well  as  to  a  poor  thermal  coupling  of 
the  generator  with  its  environment  [4].  The  output  voltage  of  the 
thermoelectric  generator  (TEG)  thus  remains  fairly  low.  In  addi¬ 
tion,  the  number  of  thermocouples  in  series  is  limited  by  the  inte¬ 
gration  constraints  and  available  area.  For  all  these  reasons,  a 
matching  converter  is  necessary  to  increase  the  voltage  as  a  means 
of  connecting  the  generator  to  a  storage  element,  in  which  voltage 
typically  lies  between  1  and  4  V  depending  on  the  technology 
introduced.  Fig.  1  displays  a  diagram  of  the  system  under  study. 

The  sizing  of  such  a  device  should  be  performed  according  to  a 
system-level  approach  that  takes  both  the  thermal  coupling  with 
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the  environment  and  the  electrical  coupling  with  the  DC-DC  con¬ 
verter  (including  losses)  into  account.  In  this  paper  however,  we 
will  simply  address  the  sizing  of  a  TEG,  potentially  with  a  heatsink, 
not  connected  to  a  DC-DC  converter.  More  specifically,  we  will  fo¬ 
cus  on  the  sizing  optimization  of  a  TEG  in  terms  of  volume  and  fol¬ 
low  this  step  up  by  presenting  a  method  that  can  be  extended  to 
other  criteria.  The  dual  thermal  and  electrical  matching  for  maxi¬ 
mized  energy  harvesting  will  also  be  studied. 

2.  Thermoelectric  system  modeling 

In  order  to  optimize  the  sizing  of  a  TEG  with  a  potential  heat¬ 
sink,  generic  sizing  models,  based  on  existing  technologies,  must 
first  be  established.  The  application  studied  herein  involves  har¬ 
vesting  heat  from  the  human  body,  though  the  models  presented 
are  easily  transposable  to  other  applications. 

2.1.  Thermal  and  electrical  models  of  a  thermoelectric  generator  (TEG) 

In  the  particular  context  of  harvesting  energy  from  the  human 
body,  the  thermal  TEG  model,  as  simplified  for  small  temperature 
differences,  is  shown  in  Fig.  2.  The  Joule  losses  and  heat  flow  due 
to  the  Seebeck  effect  can  both  be  considered  negligible,  in  compar¬ 
ison  with  the  main  heat  flow.  This  assumption  has  been  justified  in 
[4]. 

The  distinctive  features  of  such  a  system  within  the  human 
environment  are  the  very  limited  amount  of  resources  available 
and  the  poor  thermal  coupling  with  the  environment,  thus  result¬ 
ing  in  a  small  effective  difference  in  temperature  between  hot  and 
cold  heat  sources.  Given  the  level  of  recoverable  power,  the  human 
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Ambient  air:  Ta=0°C  to  30°C 


Fig.  1.  Thermoelectric  generation  system  using  the  human  body. 


Hot  source  Thermoelectric  module  Cold  source 


body  and  its  environment  can  be  considered  as  “infinite  sources  of 
temperature”,  i.e.  unaffected  by  the  presence  of  the  TEG.  Couplings 
of  the  hot  and  cold  sides  with  temperatures  Tb  and  Ta  are  repre¬ 
sented  by  two  thermal  resistances  RthB  and  RthA,  with  Rth B  being 
the  thermal  resistance  between  the  hot  heat  source  (the  human 
body  in  this  case)  and  the  hot  side  of  the  module.  The  value  of  Rt hB 
equals  the  sum  of  the  thermal  conduction  resistance  of  both  the 
human  skin  and  the  thermal  contact  between  the  skin  surface 
and  the  hot  side  of  the  module.  For  its  part,  Rt hA  is  the  thermal 
resistance  between  the  cold  side  of  the  module  and  the  surround¬ 
ing  air;  it  corresponds  to  radiation  and  natural  convection  of  the 
cold  side  with  ambient  air  and  is  very  sensitive  to  both  position 
and  movements.  Note  that  a  pin  fin  heatsink  can  also  be  added 
to  improve  heat  transfer,  but  this  would  increase  the  overall  vol¬ 
ume.  In  this  case,  we  denote  Rt hH  as  the  thermal  resistance  of  the 
heatsink.  In  addition  to  these  two  resistances,  the  thermal  resis¬ 
tance  ftthc  represents  heat  conduction  through  all  Nth  thermocou¬ 
ples,  which  are  thermally  connected  in  parallel  and  constitute 
the  TEG.  Lastly,  AT0  denotes  the  temperature  difference  between 
the  hot  source  Th  and  cold  source  Ta,  with  the  temperature  differ¬ 
ence  ATg  on  the  TEG  thermal  surfaces  being  expressed  as: 


ATg  = 


_ ftthG _ 

^thG  +  ^thB  +  ^thA 


AT0 


0) 


From  a  purely  electrical  standpoint,  the  Nth  thermocouples  are 
electrically  connected  in  series  and  therefore  display  a  total  electri¬ 
cal  resistance  denoted  Rc.  The  TEG  thus  behaves  like  an  electromo¬ 
tive  force  (emf)  Ec  (proportional  to  both  ATC  and  the  Seebeck 
coefficient  a),  associated  with  an  internal  resistance  Rc  (Fig.  3). 
Note  that  a  =  Nth  *  a0,  with  a0  representing  the  Seebeck  coefficient 
of  a  thermocouple. 

In  considering  the  notations  indicated  in  Fig.  4,  as  denoted  p 
and  X  respectively  for  the  electrical  resistivity  and  thermal  conduc¬ 
tivity  of  a  thermocouple  (average  of  the  two  materials),  hA  and  hB 
the  heat  transfer  coefficients  on  both  sides  of  the  TEG,  and  kf  the 
filling  factor  of  thermocouples  in  the  TEG,  we  can  now  express 
the  thermal  and  electrical  resistances  previously  defined  by: 


Thermoelectric  module  Electric  load 


Thermocouple 


Fig.  4.  Indications  of  the  TEG  notations. 


RthB  -  %sth ;  RthG  -  %Sth  (2) 

KthA  =  ^;  Kc  =  N^ 

Regardless  of  the  electrical  load  RL  connected  to  the  TEG,  the 
recoverable  electrical  power  Pe  can  thus  be  expressed  as  follows: 

_  kcE2c  _  kc(uATcf 

e  Rc  Rc  (  ) 

The  coefficient  kc  depends  on  the  electrical  operating  point  of 
the  TEG  and  lies  between  0  <  kc  ^  !4  For  a  perfect  electrical 
impedance  matching,  the  coefficient  value  equals  V4,  in  which  case 
the  electrical  power  is  maximized  [5,6];  this  maximum  power  is 
denoted  PeM. 
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Table  1 

Parameter  values  associated  with  TEG  and  its  environment. 


Nth 

450 

Sth 

30  cm2 

ao 

267  qV  K  1 

kh 

1.5  mm 

2 

0.77  Wm”1  K1 

Ta 

22  °C 

P 

20  qQ  m 

Tb 

37  °C 

kf 

0.6 

hA 

13  W  m-2  K-1 

kc 

0.25 

25  W  ITT2  K"1 

Fig.  5.  Evolution  of  the  maximum  recoverable  power  in  both  measured  data  (dots 
and  squares)  and  theoretical  data  (solid  and  dashed  lines)  vs.  the  number  of  stacked 
modules. 

These  benchmark  models  have  been  validated  from  experi¬ 
ments  performed  using  the  thermoelectric  module  TM-450-0.8- 
3.0  produced  by  the  Ferrotec  company.  This  module  contains 
Nth  =  450  bismuth  telluride  thermocouples,  for  a  total  surface  area 
of  approx.  Sth  =  30cm2  and  a  leg  length  of  !th  =  1.5mm.  Experi¬ 
ments  were  conducted  at  an  ambient  temperature  of  Ta  =  22  °C 
for  various  cases.  First,  the  TEG  was  placed  on  a  hot  plate  at  Th 
=  36  °C  in  considering  a  perfect  contact  between  the  TEG  and  the 
hot  plate  (fithB  =  0).  Next,  the  TEG  was  placed  in  direct  contact  with 
the  skin  of  a  human  being  (imperfect  flatness).  Moreover,  a  solu¬ 
tion  to  improve  energy  harvesting,  with  a  given  a  surface  area,  con¬ 
sists  of  stacking  identical  thermoelectric  modules  and  connecting 
them  electrically  in  series  [3,4].  For  the  current  experimental  cam¬ 
paign,  one,  two  or  three  thermoelectric  modules  were  therefore 
stacked.  The  stacking  of  thermoelectric  modules  ultimately 
amounts  to  increasing  the  leg  length  Zth  of  the  TEG,  in  the  case 
where  the  effects  of  a  module’s  alumina  plates  are  neglected  and 
only  a  one-way  heat  flow  is  considered.  It  was  also  assumed  that 
the  contact  between  two  stacked  modules  is  perfect,  which  implies 
that  the  electrical  and  thermal  contact  resistances,  as  well  as  the 
device  substrate  thermal  resistances,  can  be  neglected.  This 
assumption  was  made  due  to  the  inability  of  identifying  the  given 
resistances.  The  consequence  of  such  an  assumption  therefore  is 
that  the  thermal  resistance  of  the  TEG  RthG  is  overestimated  as 
the  leg  length  increases. 

By  taking  into  account  the  values  listed  in  Table  1,  Fig.  5  shows 
the  maximum  recoverable  power  vs.  the  number  of  stacked  mod¬ 
ules.  The  dots  and  squares  in  this  figure  correspond  to  measured 
data,  while  the  solid  and  dashed  lines  depict  the  model  dataset. 


Table  2 

Parameter  values  for  the  heatsink  model. 

km  19.5  W  m-2  K-1 

/<H2  0.67 


both  the  outer  exchange  surface  area  Sout  (i.e.  perimeter  of  the 
heatsink  multiplied  by  its  height)  and  the  inner  exchange  surface 
area  5in,  which  is  proportional  to  Sth-  In  denoting  km  and  kH2  as 
two  constant  coefficients  describing  the  heatsink  technology,  the 
model  can  now  be  expressed  as  follows: 


^Hl^th.tot  fall  (Sout  +  Sin)  /<H1  (4/ld\/Sth  + 

We  examined  the  heatsink  datasheets  provided  by  the  Aavid 
company  in  order  to  build  a  thermal  model.  The  fin  geometry  of 
such  a  heatsink  is  shown  in  the  design  diagram  in  Fig.  4.  Fig.  6  pre¬ 
sents,  for  various  fin  heights  hd,  the  variation  in  Rth H  vs.  surface 
area  Sth-  The  dots  correspond  to  values  listed  on  the  datasheet, 
while  the  solid  lines  refer  to  the  model  dataset,  with  the  values 
of  km  and  kH2  defined  in  Table  2.  (Note:  This  model  has  been 
extrapolated  since  the  max.  surface  area  of  the  heatsinks  produced 
by  Aavid  is  only  16  cm2.) 

Let  us  note  that  model  continuity  has  been  preserved.  When  the 
heatsink  height  hd  equals  zero,  the  thermal  resistance  Rth H  indeed 
equals  the  thermal  resistance  Rt^At  which  corresponds  to  radiation 
and  natural  convection,  i.e.  hA  =  km  *  kH2. 

3.  Maximizing  the  electrical  power  of  the  TEG 

According  to  the  classical  maximization  step  of  thermoelectric 
power  theory  [2],  heat  flow  through  the  TEG  is  assumed  to  be  con¬ 
stant  so  that  maximizing  TEG  efficiency  (i.e.  ratio  of  electrical 
power  to  the  heat  flow  collected  at  the  hot  source)  is  equivalent 
to  maximizing  electrical  power.  This  hypothesis  however  is  not 
applicable  to  our  study  since  heat  flow  has  not  been  imposed.  Max¬ 
imizing  the  recovered  electrical  power  [8,9]  must  then  be  treated 
as  a  different  problem  from  that  of  maximizing  efficiency,  by  use 
of  more  appropriate  means. 


2.2.  Heatsink  model 

Adding  a  heatsink  to  this  set-up  enhances  the  thermal  coupling 
between  the  cold  side  of  the  TEG  and  ambient  air  [7].  The  thermal 
resistance  of  this  heatsink  RthH  depends  on  its  total  exchange  sur¬ 
face  area  Sth_tot,  which  itself  depends  on  its  associated  capture  sur¬ 
face  area  5th  and  height  hd  (Fig.  4).  5th_tot  actually  corresponds  to 


3.1.  Thermal  impedance  matching 

The  thermal  and  electrical  models  presented  above,  in  associa¬ 
tion  with  the  heatsink  model,  enable  us  to  express  the  emf  EG  and 
electrical  power  Pe  of  the  TEG  as  follows: 

Eg  =  Ntha0AT0 ,  |th.  (5) 

kh  +  Kenv 
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— Without  heatsink  (h  =  0) 


--With  heatsink,  h  =  1  cm 
d 

--With  heatsink,  h  =  2  cm 

d 

—With  heatsink,  h,  =  5  cm 

d 

With  heat  sink,  h .  =  20  cm 

d 


30 


Leg  length  lth  (cm) 


Fig.  7.  Evolution  of  emf  EG  vs.  leg  length  /th,  both  with  and  without  a  heatsink. 


1.5 


— Without  heatsink  (h 

CL 

II 

o 

—With  heatsink,  h^  = 
a 

1  cm 

—With  heatsink,  hd  = 

2  cm 

—With  heatsink,  h .  = 

a 

5  cm 

With  heatsink,  h .  = 

d 

20  cm 

Leg  length  lth  (cm) 

Fig.  8.  Evolution  of  power  PeM  vs.  leg  length  lth,  both  with  and  without  a  heatsink. 

kck/C^ATl  Sth)th 


Pe  =  - 


(7th  +  fcei 


(6) 


The  value  of  coefficient  kenv  depends  on  the  height  of  the  heat¬ 
sink,  whereby  a  zero  value  corresponds  to  a  TEG  without  a  heat¬ 
sink.  This  coefficient  is  expressed  below  in  taking  into  account 
the  previous  thermal  models: 


^env  — 


1 


\ 


By  inputting  the  parameter  values  listed  in  Tables  1  and  2,  Figs. 
7  and  8  show  respectively  the  evolution  of  emf  Ec  and  electrical 
power  PeM  on  the  adapted  load  (/<c  =  %),  both  with  and  without  a 
heatsink  (i.e.  hd  =  0),  as  a  function  of  the  leg  length  lth. 

Given  a  surface  area  Sth  and  leg  length  Zth,  the  addition  of  a  heat¬ 
sink  can  enhance  the  emf  EG  as  well  as  the  electrical  power  PeM.  The 
shape  of  PeM  is  known  [10]  and  exhibits  an  optimum  value  for  leg 
length  as  a  ratio  of  heatsink  height,  though  this  optimal  length 
tends  to  be  technologically  infeasible  (i.e.  too  high).  One  approach 
to  achieving  basically  the  same  improvement  would  be  to  stack  the 
modules,  as  pointed  out  above. 

A  fundamental  problem  of  harvesting  energy  from  the  human 
body  relates  to  the  poor  thermal  matching  of  a  TEG  with  its  envi¬ 
ronment.  In  other  words,  the  thermal  resistance  Rth G  is  generally 
very  low  compared  to  the  environment-related  resistances  PthB 
and  PthA  (or  RthH  with  a  heatsink).  Let  us  denote  RthE  as  the  sum 
of  thermal  resistances  related  to  the  environment,  i.e.  RthB  +  PthA 
without  a  heatsink  and  RthB  +  RthH  with  a  heatsink.  Since  the  inter¬ 
nal  thermal  and  electrical  resistances  of  the  TEG,  denoted  respec¬ 
tively  flthc  and  RG,  are  proportional  to  the  leg  length  Zth,  an 
increase  in  leg  length  serves  to  improve  the  thermal  coupling 


and  thus  increases  ATG,  yet  also  enhances  the  internal  resistance 
Rg  and  therefore  a  priori  decreases  the  maximum  electrical  power 
harvested.  In  considering  Eq.  (6),  it  is  possible  to  prove  that  electri¬ 
cal  power  is  maximized  for  an  optimal  value  of  Zth  [11],  which  de¬ 
pends  on  heatsink  height  and  thus  differs  whether  or  not  the 
heatsink  is  present.  This  condition  is  expressed  as  follows: 


/thopt  —  Akf  |  ^ 


1 


\ 


km 


(8) 


The  maximum  power  is  obtained  when  the  internal  thermal 
resistance  of  TEG  RthG  equals  the  sum  of  thermal  resistances  in 
its  environment,  i.e.  summing  to  RthE.  The  difference  in  tempera¬ 
ture  A Tg  across  the  TEG  is  thus  equal  to  half  the  difference  in  tem¬ 
perature  between  the  body  and  ambient  air  [12],  a  situation  that 
corresponds  to  thermal  impedance  matching. 

Increasing  the  leg  length  (or  stacking  thermoelectric  modules) 
thus  improves  thermal  coupling  for  a  given  surface  area,  which  in 
turn  increases  the  emf  and  electrical  power  of  the  TEG.  This  result 
proves  to  be  a  considerable  advantage  when  connecting  the  TEG 
to  a  DC-DC  converter.  At  these  power  levels,  major  losses  in  the 
converter  are  in  fact  switching  losses  that  decrease  as  input  voltage 
increases  [13].  The  best  solution  to  reduce  these  losses  would  be  to 
include  a  high  number  of  thermocouples  in  the  TEG;  for  example, 
the  TEG  produced  by  the  Micropelt  and  Thermolife  companies  are 
made  of  a  thin  film.  This  technology  allows  integrating  a  greater 
thermocouple  density  (approx.  5000  per  cm2  compared  to  15  per 
cm2  for  the  module  TM-450-0.8-3.0  produced  by  the  Ferrotec  com¬ 
pany)  and  therefore  yields  high  open  circuit  voltages  EG  of  TEGs  rel¬ 
ative  to  the  surface  area  of  the  hot/cold  plate.  The  solution 
discussed  above  of  stacking  modules  also  offers  a  way  to  increase 
Eg  when  the  leg  length  is  limited  for  technological  reasons. 

Increasing  the  leg  length  and/or  adding  a  heatsink  however  also 
expands  the  size,  i.e.  volume,  of  the  TEG,  which  can  now  be  ex¬ 
pressed  as: 


^th  —  ^th  (/th  +  frd) 


(9) 


As  a  result,  it  may  be  worthwhile  to  size  the  TEG  by  seeking  to 
minimize  volume  while  maximizing  energy  recovery;  this  ap¬ 
proach  will  be  discussed  further  below. 


(7)  3.2.  System  sizing  optimization 


In  order  to  size  a  TEG  that  maximizes  the  electrical  power  pro¬ 
duced  while  minimizing  its  volume,  we  have  introduced  a  bi-objec- 
tive  multi-variable  optimization  algorithm  based  on  particle  swarm 
optimization  [14]  and  Pareto  dominance.  This  algorithm,  called 
MOPSO,  has  been  implemented  in  a  Matlab  environment  in  our  lab¬ 
oratory  and  has  already  been  used  for  another  application  [15].  We 
were  thus  able  to  optimize  TEG  sizing,  both  with  and  without  a 
heatsink,  in  considering  electrical  impedance  matching  (i.e.  the 
TEG  voltage  is  thus  equal  to  half  its  emf).  With  a  given  number  of 
optimization  parameters  (here  5th,  lth  and  hd  if  a  heatsink  is  pres¬ 
ent),  this  optimization  algorithm  returned,  after  several  iterations, 
a  set  of  optimal  solutions  meeting  the  input  specifications  and  dis¬ 
tributed  on  a  Pareto  front.  The  two  competing  criteria  chosen  are: 
minimization  of  the  TEG  and  heatsink  (if  present)  volume  Vth;  and 
maximization  of  harvested  electrical  power  PeM.  In  addition,  we 
forced  the  geometric  parameters  to  vary  only  within  a  bounded 
interval  (see  Table  3).  Note  that  the  constraint  introduced  to  limit 
heatsink  height  to  10  cm  is  arbitrary,  but  it  was  necessary  to  set  a 
limit  in  order  to  respect  the  validity  of  this  thermal  resistance  mod¬ 
el.  For  these  optimizations,  the  temperature  values  of  the  body  ( Tb ) 
and  ambient  air  (Ta)  are  those  listed  in  Table  1 .  Moreover,  a  constant 
thermocouple  density,  equal  to  that  of  module  TM-450-0.8-3.0 
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Table  3 

Bounded  interval  of  the  geometric  parameters  for  the 
optimization  algorithm. 


Parameter 

Min 

Max 

Sth 

0 

100  cm2 

Zth 

0 

10  cm 

hd 

0 

10  cm 

(Ferrotec),  was  assumed  due  to  a  lack  of  sufficient  generic  informa¬ 
tion  on  the  effects  of  thermocouple  density.  The  number  of  thermo¬ 
couples  Nth  is  thus  proportional  to  the  surface  area  Sth.  Lastly,  let  us 
recall  that  it  was  decided  to  vary  the  leg  length,  which  is  roughly 
equivalent  to  stacking  modules  (with  an  exception  made  for  the 
thickness  of  the  modules’  alumina  plates). 

Fig.  9  shows  the  Pareto  front  solutions  of  optimizations  with 
volume  minimization  and  power  maximization.  The  red  dots  cor¬ 
respond  to  the  optimization  of  a  TEG  with  a  heatsink,  while  the 
blue  dots  reflect  the  TEG  optimization  without  a  heatsink.  For  this 
second  optimization,  we  forced  the  value  of  parameter  hd ,  which 
corresponds  to  heatsink  height,  to  zero,  i.e.  hd  =  0. 

It  can  be  noticed  that  the  two  Pareto  fronts  obtained  overlap  at 
first.  Then,  beyond  the  coordinate  point  Vth  =  120cm3  and 
PeM  =  0.8  mW,  the  effect  of  adding  a  heatsink,  given  a  certain  vol¬ 
ume,  makes  it  possible  to  harvest  greater  power. 

As  a  preferred  approach  to  interpreting  these  results,  we  stud¬ 
ied  the  evolution  in  optimization  parameters,  i.e.  surface  area  Sth 
(Fig.  10),  leg  length  lth  and  heatsink  height  hd  (Fig.  11),  as  well  as 
the  evolution  in  thermal  resistances  of  both  the  generator  Rth G 
and  the  environment  PthE  (Fig.  12). 

For  the  sizing  optimization  of  a  TEG  that  harvests  maximum 
power  as  a  result  of  electrical  impedance  matching,  let  us  observe 


Fig.  11.  Evolution  in  the  optimal  leg  length  lth  and  heatsink  height  hd  vs.  electrical 
power  PeM- 


Electrical  power  PeM  (W)  x  1 0 


Fig.  12.  Evolution  in  thermal  resistances,  RthG  and  PthE,  vs.  electrical  power  PeM. 


Fig.  9.  Bi-objective  optimization  solutions  that  minimize  volume  and  maximize 
productivity. 
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Fig.  10.  Evolution  in  the  optimal  surface  area  Sth  vs.  electrical  power  PeM- 


that  without  the  presence  of  a  heatsink,  the  surface  area  5th  chosen 
by  the  optimization  algorithm  is  always  the  maximum  allowed,  i.e. 
here  100  cm2  (blue  dots  in  Fig.  10).  According  to  Eqs.  (6)  and  (7), 
regardless  of  the  electrical  operating  point  of  the  TEG,  its  power 
is  always  proportional  to  Sth-  Hence,  in  order  to  increase  the  har¬ 
vested  power  with  a  given  surface  area,  the  leg  length  Zth  must  sim¬ 
ply  be  increased  to  its  optimum  value  (see  Eq.  (8)).  In  Fig.  11,  the 
evolution  of  Zth  without  a  heatsink  basically  corresponds  to  the 
increasing  part  of  the  PeM  vs.  Zt h  graph,  as  plotted  in  blue1  in 
Fig.  8.  Once  the  maximum  power  has  been  reached,  the  increase 
in  Zth  is  no  longer  beneficial.  For  a  desired  electrical  power  and  gi¬ 
ven  heat  exchange  conditions,  it  then  becomes  preferable  to  intro¬ 
duce  large  surface  areas  and  a  thin  TEG.  Moreover,  in  Fig.  12,  the 
evolution  in  thermal  resistances  of  both  the  generator  Rth G  and 
environment  Rth E  confirms  that  maximum  power  is  achieved  for 
a  satisfactory  thermal  impedance  matching,  i.e.  when  Rth G  =  RthE- 
For  the  sizing  optimization  of  a  TEG  including  a  heatsink,  the 
electrical  power  is  no  longer  proportional  to  surface  area  Sth  given 
that  coefficient  kenv  depends  on  Sth ,  as  a  result  of  the  heatsink  mod¬ 
el  (Eq.  (7)).  The  value  of  Sth  selected  by  the  algorithm  (red  dots  in 
Fig.  10)  however  remains  the  maximum  allowed.  Regarding  the 
evolution  of  Zth  and  hd ,  Fig.  11  shows  that  at  low  power  levels, 
the  heatsink  height  hd  equals  zero,  which  means  that  adding  a 
heatsink  is  not  advantageous.  Thanks  to  the  continuity  of  the  heat¬ 
sink  model  (Eq.  (4)),  the  two  Pareto  fronts  illustrated  in  Fig.  9  can 
thus  be  combined.  Fig.  1 1  also  shows  that  at  higher  power,  the  TEG 
with  a  heatsink  height  greater  than  the  leg  length  is  to  be  pre¬ 
ferred.  In  addition,  it  has  been  observed  that  the  power  attained 


1  Please  note  that  Figs.  8  and  10  will  appear  in  B/W  in  print  and  color  in  the  web 

version.  Based  on  this,  please  approve  the  footnote  1  which  explains  this. 
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is  higher  when  adding  a  heatsink,  since  power  increases  with  heat¬ 
sink  height  hd  (Eqs.  (6)  and  (7)).  Maximum  power  is  actually  at¬ 
tained  by  the  algorithm  when  the  heatsink  height  equals  the 
maximum  allowed  value,  i.e.  10  cm  (Fig.  11).  Lastly,  thermal 
impedance  matching  has  been  observed  for  a  TEG  with  a  heatsink 
when  electrical  power  is  maximized  (Fig.  12). 

In  Figs.  13  and  14,  the  power  density  (i.e.  power  per  unit  vol¬ 
ume)  has  been  plotted  vs.  both  electrical  power  and  volume, 
respectively. 

Let  us  note  in  Figs.  13  and  14  that  whether  with  or  without  a 
heatsink,  volume  increases  faster  than  electrical  power  (i.e.  the 
power  density  function  is  decreasing). 

Furthermore,  these  optimizations  could  be  conducted  by  setting 
power  density  maximization  as  the  optimization  criterion  rather 
than  electrical  power  maximization.  We  opted  for  this  approach 
and  found  the  same  results  as  presented  above. 

4.  Conclusion 

In  this  paper,  we  have  presented  the  sizing  optimization,  in 
terms  of  volume,  of  a  thermoelectric  generator,  which  as  an  option 
may  include  a  heatsink,  in  the  specific  context  of  harvesting  energy 
from  the  human  body. 

Models  of  a  thermometric  module,  simplified  for  small  temper¬ 
ature  differences,  and  a  heatsink  were  built  from  experimental  val¬ 
idation  or  datasheets,  and  both  types  of  models  were  based  on 
existing  technologies.  Moreover,  continuity  between  the  two  mod¬ 
els  has  been  well  respected,  meaning  that  for  a  given  surface  area 
the  thermal  resistance  of  a  heatsink  with  zero  height  corresponds 
to  the  radiation  and  natural  convection  resistance  of  the  thermo¬ 
electric  module.  This  study  of  TEG,  with  or  without  a  heatsink, 
highlights  the  relevance  of  electrical  and  thermal  matching  of  the 
TEG  to  its  environment  in  order  to  maximize  its  power.  Thermal 


impedance  matching  can  be  performed  by  increasing  the  TEG  leg 
length  until  its  optimum  value  or,  should  this  prove  technologi¬ 
cally  infeasible,  by  stacking  thermoelectric  modules. 

The  models  developed  models  have  enabled  the  sizing  optimi¬ 
zation  of  such  a  thermoelectric  generator.  By  defining  two  compet¬ 
ing  optimization  criteria,  consisting  of  minimizing  the  TEG  volume 
(with  the  inclusion  of  a  potential  heatsink)  while  maximizing  the 
power  produced  by  the  generator,  we  have  demonstrated  that 
from  a  certain  volume,  adding  a  heatsink  makes  it  possible  to  har¬ 
vest  more  power.  Moreover,  we  have  shown  that  the  solution  min¬ 
imizing  volume  leads  to  an  ultra-thin  TEG  with  power  being 
maximized  when  TEG  impedance  is  matched  thermally  and  elec¬ 
trically  with  its  environment.  If  the  capture  surface  area  is  indeed 
limited,  then  the  minimization  of  total  volume  yields  a  greater 
thermocouple  leg  length  (e.g.  by  stacking  modules).  Under  these 
conditions,  maximum  electrical  power  is  obtained  for  a  leg  length 
corresponding  to  thermal  matching  (equality  between  the  TEG 
thermal  resistance  and  the  thermal  resistance  related  to  the  envi¬ 
ronment).  At  higher  electrical  power,  a  heatsink  is  necessary  to  im¬ 
prove  thermal  coupling  with  the  environment,  and  the  total 
volume  is  still  minimized  for  a  leg  length  (smaller  this  time)  lead¬ 
ing  to  thermal  matching  with  a  heatsink. 

Furthermore,  the  TEG  studied  herein  displays  a  low  thermocou¬ 
ple  density.  Technologies  using  higher  density  (though  a  priori  the¬ 
oretically  less  productive  per  unit  area)  could  be  better  suited  to 
such  a  system  optimization.  It  might  prove  beneficial  to  define 
thermocouple  density  as  an  optimization  parameter. 

To  the  best  of  our  knowledge,  such  global  optimization  results 
are  entirely  new,  with  this  method  being  potentially  extended  to 
other  criteria  as  well,  i.e.  economic  or  environmental,  provided  a 
sufficient  amount  of  data  has  been  associated  with  these  criteria. 
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